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Abstract. Modeling plasma physical processes in astrophysical context demands for both detailed 
kinetics and large scale development of the electromagnetic field densities. 

We present a new framework for modeling plasma physics of hot tenuous plasmas by a two-split 
scheme, in which the large scale fields are modeled by means of a particle-in-cell (PIC) code, and 
in which binary collision processes and single-particle processes are modeled through a Monte- 
Carlo approach. Our novel simulation tool - the PhotonPlasma code - is a unique hybrid model; 
it combines a highly parallelized (Vlasov) particle-in-cell approach with continuous weighting of 
particles and a sub-Debye Monte-Carlo binary particle interaction framework. 

As an illustration of the capabilities we present results from a numerical study [ 1 ] of gamma- 
ray burst - circumburst medium interaction and plasma preconditioning via Compton scattering. 
We argue that important microphysical processes can only viably be investigated by means of such 
"trans-Debye scale" hybrid codes. 

Our first results from 3D simulations with this new simulation tool suggest that magnetic fields 
and plasma filaments are created in the wakefield of prompt gamma-ray bursts. Furthermore, the 
photon flux density gradient impacts on particle acceleration in the burst head and wakefield. 
We discuss some possible implications of the circumburst medium being preconditioned for a 
trailing afterglow shock front. We also discuss important improvements for future studies of GRB 
wakefields processes, using the PhotonPlasma code. 
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1. INTRODUCTION 

The prompt emission (PE) from a gamma-ray burst (GRB) precedes the GRB ejecta 
from the central engine by seconds-to-hours (depending on the phase in consideration). 
The PE will strongly influence the circumburst medium (CBM) before the arrival of the 
relativistically trailing shocked ejecta that subsequently sweeps up and shocks the CBM. 
We have sketched this scenario in figured] 

Lorentz factors of the GRB ejecta are generally theorized to be much higher (of order 
r « 10 2 - 10 3 , e.g. y, |j, |4fl) than that for any motion of the GRB progenitor - and a 
possible accompanying wind - itself. Both interpretation of observations and modeling 
of relativistic shock physics and emission, concerned with GRB early-to-late afterglows, 
is complicated by its relativistic nature. The plasma microphysics of GRB afterglow 
shocks propagating through a CBM has been studied and extensively modeled using 
particle-in-cell (PIC) models in recent years [ULSEzLlsLl^l]. 

In contrast, modeling and observational interpretation of the PE-CBM interaction 



is somewhat simplified; whereas the GRB ejecta propagates at some (unknown) high 
Lorentz factor and involve a strong dependence on initial burst conditions, rather, the 
CBM is likely to be at rest relative to the progenitor from whence the burst emanates. 
To good approximation, this is true for the entire duration of the prompt flash. CBM 
rest frame events and quantities will - after cosmological corrections - be directly 
translatable to the BATSE rest frame, and since the PE is the first to reach the observer 
it is a strong indicator on what environment the trailing ejecta will experience. We 
can more directly interpret observational data and PE-CBM interaction modeling, even 
while relying on fewer assumptions. Unknowns are reduced to a few: CBM density, no, 
prompt emission flux through a CBM surfaceQ, and the injection photon spectrum from 
the prompt fireball emission (v,F v (t)). 
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FIGURE 1. (Color online) Rendering - and not to scale: Prompt flash emission from a GRB central 
engine precedes a shocked ejecta, and traverses (moving right) a quiescent CBM. After the prompt 
gamma-ray front has passed, the CBM will relax until the ejecta shock front sweeps up the preconditioned 
CBM plasma. Clearly, the prompt emission front continues all the way back to the ejecta, although photon 
energies decrease with decreasing distance from the shocked ejecta. 

To investigate stochastically induced wakefield effects, one may in a general way pose 
the following central questions concerning the PE-CBM interaction: 

• How does the PE pulse influence the CBM plasma. In particular; how does it affect 
particle populations and electromagnetic fields? 

• How does the CBM plasma affect the PE gamma-ray pulse? 

• Which physical processes are needed to adequately describe a prompt y-ray burst 
signal as it traverses a circumburst plasma? 

The rapid variability in the prompt signal has been a matter of hefty debate since the 
advent of GRB research in the late 1960's. The community has been arguing in favor 
internal shock models on one hand, and external shocks models on the other. The same 
variability has proven to be a cornerstone in establishing the relativistic expansion aspect 
ofGRBs. 

Modeling the relevant micro-physics and plasma effects is central to an understanding 
of whether the variability observed in the prompt phase (e.g. by BATSE) and - more 
recently - the variability observed in the early afterglow phase (SWIFT), should be 
attributed to an internal or external shock model, and perhaps if we need to revise these 
models altogether. 



We simplify by assuming a slab geometry for the microphysical scenario 



2. MOTIVATING A CHOICE OF MODEL 



We model, for the first time |Q]], the detailed interaction between a Vlasov plasma (see 
discussion below, and Eqn[T|) and a photon burst through stochastic wakefield processes, 
motivated by the scenario of a prompt gamma-ray burst propagating through a quiescent 
circumburst medium plasma [10]. For this purpose, we are using a unique and highly 
improved particle-in-cell scheme that contains a framework to handle both Vlasov 
plasmas through a particle-in-cell (PIC) code described elsewhere [|llL LL20, and sub- 
Debye sphere binary particle processes. For modeling the GRB-CBM interaction we 
employ only Compton scattering - which is also described elsewhere ||l2Lll3|l . 



2.1. Trans-Debye Kinetic Plasma Modeling 

Photons with energies in the BATSE energy window, 10 4 eV < E r < 10 6 eV, in- 
volved in scattering processes during the very first GRB-CBM interaction, will have 
wavelengths much smaller than: 1) the plasma interparticle distance and 2) any charac- 
teristic plasma phenomenon. Consequently, any encounter between typical plasma par- 
ticles (electrons or ions) and photons can be viewed as a binary collision taking place on 
time and spatial scales well below the characteristic plasma scales. 

Furthermore, the two plasma parameters central to a kinetic description of plasmas^ 
are the Debye length, Ad, and the electron skin-depth, 8 e . Their ratio is given by 
Ap I 8 e ~ 1.4- 10~ 3 r 1 / 2 , where the temperature, T, is expressed in eV uM- This 
ratio is less than one for all reasonable CBM temperatures, and tends to one only for 
very high electron temperatures, T > 10 6 eV. For the circumburst medium of GRB, it is 
much less than unity; for the post-flash CBM the ratio may approach unity. 

Conventional Particle-In-Cell codes. Fundamentally PIC codes implement the 
Vlasov approximation; particles interact only through their self-consistently generated 
smoothed fields, and collisions are only of a wave-particle nature. Adding collision 
effects - other than wave-single-particle effects - ultimately leads to the Vlasov- 
Boltzmann formulation for plasmas 

dtf + V ■ d r f+ -(E + v x B) • d p f = [d t f] coll . Vl =° V , (1) 
m s 

where the electromagnetic forces/sources are obtained by interpolating to/from particles 
to the computational grid that defines the electromagnetic field. The fields are in turn 
integrated using Maxwell's equations, and particles are moved under the (relativistic) 
Lorentz force. An lower limit may be imposed on Ap; although PIC codes are ideally 
collisionless, integration errors will introduce artificial collisions in the Vlasov-plasma, 
[dtf]coll. 7^ 0; the plasma heats (from phase space diffusion) until the Debye length is 



slightly larger than the grid cell length II 1511 ; Ad 1.5Az. For other reasons, concerning 



Un-magnetized or weakly magnetized plasmas. 



artificial damping of collectively exicted phenomena iflll section 3.2ff], the Debye 
length should be chosen not too far away from the grid scale; a stricter limit on Ad 
results. Consequently, we are constrained by generic particle-in-cell code properties to 
choose Ax: ~ Ad and all particle-particle processes such as scattering, pair production 
and decay are as argued all defined well below the grid scale. 

The Random Phase Approximation. For a plasma consisting of an ensemble of point 
particles, p(r) = £,-5(r,- — ro), an equation of motion (e.o.m.) can be written for the 
electron density. Following Pines and Bohm [16], we write 

p,= -£(k-v ; )V'' k ' r <- £ [4xe 2 /m(k') 2 }k-k' {exp[i(k' -k) -r^expi-ik' , 

i k'ij,k'^0 

(2) 

where the e.o.m. is now written in terms of Fourier components. The last term can then 
be split into two parts, corresponding to sums over k ; = k and k' 7^ k, wherein, the first 
part reduces to the total number of particles in the ensemble. The second part - for all 
k' 7^ k - reduces to = 1 (rather, the exponential i(k ; — k) ■ r, = 0) under the assumption 
that a very large number of particles are distributed at random, spatially, and that their 
phases consequently average to zero. Pines and Bohm refer to this as the random phase 
approximation jl6\ . This approximation leads to a reduced expression for the e.o.m. for 
the density: 

p k = -£(k. Vi )V^ r ') - (W/m)£e-( k - r ') , (3) 

i i 

which for small k becomes the well known expression for electron plasma oscillations. 
For all k both terms in equation |3] must be taken into account and the relative strength 
of the two determine whether collective or single-particle effects are more important. 
Consequently, the random phase approximation ensures that for k 2 <C A^ 2 a collective 
description is suitable, whereas for k 2 ^> A D 2 the dynamics is dominated by single- 
particle effects. 

It is this cross-over or split dynamics approximation that motivates our use of the hybrid 
scheme devised in the PhotonPlasma code; the PIC code handles large-scale (large 
compared with the Debye length) plasma effects, while the Monte-Carlo collisional 
framework handles collisional effects at short wavelengths. 

2.2. Detailed Particle-Particle Processes 

Every macro-particle in PIC codes represents a large number of real physical particles, 
sometimes as many as 10 30 . Compared with previously employed PIC schemes, the 
PhotonPlasma code has been improved to include variable weighting of every single 
macro-particle. The charge density defined in terms of an interpolation of point particles 
to the field grid, for example, is then modified: 



Pc(r,0=£9oJ^5(ri-ro,0 — ► p c (r,t) = J^aJ^w/Sfo-ro,*) , (4) 

a i a i 



with Wi the weight of particles, and with q a the species' charges. This seems relatively 
benign at first, since the gain in resolution for the macroscopic quantities is real valued 
but only discretely so. The advantage becomes clear when considering scattering pro- 
cesses such as Compton scattering, p e +p y — > p e +Py, used in the GRB wakefield 
experiments presented in the following sections. With variable weighting, the physical 
number densities of the plasma constituents are rather now given in some volume, e.g. a 
cell volume, by 

n e = J>W , n Y = 2>J», (5) 

i=l j=l 

for electrons and photons respectively. We can now define the physical Compton scat- 
tering rate in that same volume as 

with Oc the Compton scattering cross section. We have assumed here that the num- 
ber of scatterings is proportional to the product number densities of the interacting 
species. This leads to a change in number density of the original interacting particles, 
An e , r oc n e ttyOc(p e ,p.y)cA£, where resolution gain is now real valued and continuous. 
The number of (Compton) scatterings is then given as a continuous real number quan- 
tity, An e j on the cell volume, even if the particles are discrete entities. 

Scattering is achieved by splitting the macro-particle pair weights into a total of 
four new particles, two of which carry the new scattered momenta and energies as 
determined by the microscopic physical cross section (CTclP^Py))- This is manifested 

through = Wnew + w scatt > f° r m e electron and photons both. An important fact 
about this scattering procedure is that we retain full resolution in change of momenta, 
energies, weights, charges decay rates etc. For example, consider the scattering in a 
volume of 10 30 photons and electrons. Effectively, for a scattering fraction of 10 , 
10 24 physical particles are scattered away from a 'mother' particle of size 10 30 . This 
illustrates the flexibility and the fundamentally continuous nature of the trans-Debye 
kinetic approach. 

For the modeling of GRB prompt emission interacting with the tenuous and relatively 
cold circumburst medium, the considerations summarized in this section motivated the 
development and deployment of the PhotonPlasma code. 

3. SETUP - STOCHASTIC WAKEFIELD SIMULATION 

A synthetic thermal GRB is delivered to the quiescent CBM plasma by adding photons 
(particles) to the computational domain on the left volume boundary, according to a 
prescribed light curve and spectrum - this is sketched in figure [2] 



Gamma-Ray Burst Synthesis. The light curve for the synthetic gamma-ray burst is 
inspired by observational data fitted to a FRED function ||17Lll8|l with parameters chosen 



for BATSE trigger 3891 (GRB951 102), r = 1.26 and d = 2.67. We have set the duration 
of the photon burst to Tqrb = 200ft>~ , which corresponds to about 1% of T$q for trigger 
3891. Assuming a pure blackbody spectrum, the accompanying temperature curve is 
shown in the lower left panel of figure [3] 



FIGURE 2. (Color online) Rendering of the prompt GRB photons (colored particles) traversing the 
quiescent CBM (grey particles) in our model. Photons are added to the computational domain on the left 
boundary. As the photon number flux density increases, the weights (size of particles) - rather than the 
number - of the photons vary in a continuous way. Color designates that photons further carry different 
energies (through v) as well. This detailed balance approach saves computational effort and increases 
speed. The microphysics is unchanged (section|2|. 

A considerable sample from the BATSE catalogue may be modeled with better statis- 
tic assuming that the prompt emission spectrum is composed of a power law component, 
overlayed with a blackbody spectral component (see e.g. lfl9ll ). Temperatures are found 
to be of order Tbb ~ 50keV. This is well below pair threshold and - in the first ap- 
proximation - we may neglect pair production as an important effect during late time 
evolution of the prompt burslH. Photon energies for the thermal evolution are normalized 
to the BATSE observational window, Ey £ [33. 6ke V ,3.3 6Me V] which is corrected for a 



FIGURE 3. (Color online) Left: synthetic GRB fast-rise-exponential-decay (FRED) light curve (top), 
and the temperature evolution (bottom) of the PE, assuming perfect blackbody radiation. Time, flux 
and temperature are normalized to the total burst duration, peak flux and temperature at peak flux - 
respectively. Right: the spectral photon radiance of the thermal burst for various times. Photon energy and 
number density is scaled in per cent of BATSE energy window and peak radiance - respectively. Color of 
individual curves designates linear time from burst onset. 

Temporal evolution of the GRB is achieved through the assignment of variable 
weighting to a constant number flux of computational macro-particle photons (see also 
section O, this corresponding to a varying number flux of physical photons, in accor- 



3 Note, however, that Svensson ll20ll found that even for these moderate temperatures, pair production 
could be significant in the high-energy tail of the particle distribution - see also 12111 . 





Synthetic Burst Evolution 




dance with the prescribed light curve (particle weights) and spectrum (particle energies). 
Photon energies are sampled from MC integration of the (time-dependent) Planckian; 
we obtain a comparatively high accuracy at low computational cost. The photon spectral 
radiance (spectral photon number density) is given by 

F(v t) Nyiy,t) Wy(v) 

« y (V,?)oc ^oc £ Wi(V,t) => n y (V,t)°c £ Wi(t) . 

V !=1 (=1 

Through split assignment we impart time- variability to the weights, Wi(v,t) — > Wi(t), 
and frequency variability to individual particles, Ny(v,t) — > Ny(v). This yields a high 
degree of flexibility compared to conventional PIC codes as explained in detail in 
section [21 The macro-particle photon number density is set constant at 40 particles/cell 
for the entire duration of the burst. 

Circumburst Medium Plasma Model and Boundary Conditions. We assume the 
CBM to be a quasi-neutral thermal hydrogen plasma (ions and electrons) of constant 
density with moderate temperature, v t h, e ~ 0.1c, and slightly reduced mass ratio of the 
plasma constituents of m,/m e = 256. The plasma is initially thermal and field free (ex- 
cept for thermal fluctuations). No kinetic photons are present in the computational vol- 
ume. Boundaries are periodic in the (X-Y-) plane transverse to the GRB pulse propaga- 
tion. A grid resolution of 100x20x4000 translates to {L x ,Ly,L z } ps {13<5 e x2<5 e x500<5 e } 
for our choices of length units and plasma density. Initially, the number of plasma par- 
ticles is 40/species/cell, but the total number of particles increases as scattering splits 
the interacting particle pairs during the run. Using a relatively flat volume aspect ratio 
is a compromise to ensure skin-depth resolution at the lowest possible grid size to give 
any 3D structure a chance of growing marginally without spending excess computing 
time. The moderate number of particles/cell has minor impact on resolution issues as 
discussed in section [2l since scattering adds particles and conserves the momentum lo- 
cally everywhere, and globally conserves energy to better than 1%. The burst duration 
in light travel length is Lqrb ~ 64 5 e . The lower/upper Z-boundary (GRB influx/outflux 
boundary) specularly reflects/absorbs all plasma, but both ends allow for photons to exit. 
Electromagnetic fields are damped in a thin layer, AL Z ~ 0. 1 %L Z on both upper and lower 
Z-boundaries and are periodic in the tranverse direction (XY), to avoid any spurious EM 
wave growth. 

4. RESULTS - STOCHASTIC GRB WAKEFIELD EFFECTS 

Here we present results only from the largest run in our wakefield simulation batch, 
with parameters as described in the previous section. For a more complete treatment 
of various effects of changes to the setup, please refer to [flo|l and HI] where also 
a number of tests and details are given. Robustness of the central results presented 
here was tested through a consistent scan of parameter space. We concentrate here on 
the plasma-particle heating/acceleration aspects and electromagnetic field growth. We 
shall distinguish between 'plasma particles' (electron macro-particles) and 'photons' 
(photonic macro-particles), see also section [2] Figures 0] and \5\ compactly summarizes 



the results of stochastically forcing the quiescent thermal plasma - the CBM - by a 
strong thermal photonic pulse - the GRB. 
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FIGURE 4. (Color online) Left - electron scatter plot: Color designates weight in log-scale - blue 
particles have weights Wblue ~ 10~ 5 wo, where wq is the initial weight. The yellow line coincides with zero 
bulk z-momentum, p z y = 0. For clarity photons (black) are offset vertically by p z< y = +6. Right - photon 
scatter plot: The initial burst of photons (in color - moving right) is visible from t ~ 1500 to t ~ 1300. 
Pink arrows: local photon upscattering. Green lines: energy threshold for y+ y ^ e + +e~ for photons 
pairwise oppositely outside (away from p z = 0) these lines. Electrons (black) are offset by p ze = +6. 
Here, p z> y = [hv] z /m e c 2 . 

Wakefield - Particle Effects. As the photon pulse interacts with the plasma through 
Compton scattering, it excites waves in the plasma. Although excited by single macro- 
particle scatterings in a highly optically thin system (t -C 1; the integrated weight of 
photons scattered at least once per wavelength is less than 10~ 3 ), these waves are the 
collective reaction of the plasma to the photon stochastic forcing. We emphasize the 
random nature of the scattering. 

A strong acceleration of the particles in the wakefield of the photon pulse is seen in 
Figure |4] - right panel, where plasma particles are plotted in color. To sustain quasi- 
neutrality in the wakefield plasma, a subset of the plasma particles (taken from those 
not scattered) are accelerated backwards, p z < 0, in the wakefield. Their response is 
collective as no preference is given to any one particle in the Vlasov-formulation; all 
processes are described by characteristic scales, kc, above the Debye-length, kc > A^ 1 . 

Plasma structures excited by the photon pulse are almost co-moving with the pulse, 
lagging gradually behind the pulse front due to differences in propagation speed of the 
photons and Langmuir waves. The structures are therefore dragged out behind the first 
photon-plasma encounter, all while growing in strength (height), and they increasingly 
interact with the photon pulse; a feedback loop to the photon pulse produces the photon 
energy spikes shown in the right panel of Figure HI 

Wakefield - Electromagnetic Field Effects. A counter- streaming plasma results from 
the anisotropic nature of the Compton interaction as seen in the left panel of Figure H 
where the momentum is centered on |pj = 0, and we may expect the development of rel- 
ativistic (or mildly so) streaming instabilities, such as nTamentation and the two-stream 
instability (see e.g. [|23l . l24ll N ). Figure[5]displays a cross section through the computational 
domain, revealing the production of current filaments (nTamentation) and an associated 



magnetic field, with an energy density equal to a few per cent of equipartition, £g « 5%, 
with the total kinetic energy transfer fed to the plasma. While ions play the role of main- 
taining the filaments, the electrons produce the current necessary for the magnetic field 
to exist. We have tested the robustness of the field generation for various burst durations 
and luminosities, and found B-field generation in all cases II 1011 . 
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FIGURE 5. (Color online) Current structures for electrons (panel 1 ) and ions (panel 3), the photon pulse 
(panel 2), and (transverse) magnetic field (panel 4) in our 3D wakefield simulation, for t = 1020<B ; 7 f , 1 w 
10 4 Af. The GRB photon pulse has light-length L c = c/x<jrb = 200 ~ 64-8 e . The computational volume 
has dimensions {L x ,L y ,L z } = {40, 8, 1600} or, in electron skin-depths, {L x ,L y ,L z } w { 125 e , 28 e , 5005 e }. 



A Photon Anisotropy Effect. The energy ex- 
change in Compton scattering is 
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Forwardly scattered photons exchange rather 
small amounts of energy with the (rest) elec- 
tron whereas reversely scattered photons may ex- 
change large amounts of energy. The figure inset 
shows differential Compton cross section Q25[| for 
energies spanning the BATSE observational en- 
ergy window. High energy photons (electron rest 
frame) will scatter mainly forward and low energy 
photons mainly either backwards or forwards. 

We may now explain the apparent stability of 
the CBM counter- streaming plasma set up by the 

GRB pulse: photon scattering on electrons moving forward (pulse-parallel) will scatter 
severely and impart energy to the electron so, as to accelerate it even more (CBM rest 
frame). For a photon scattering on an electron moving backward (pulse-anti-parallel), 
however, the momentum and energy exchange in the scattering would be such that no 
significant acceleration of the electron would result. In a second scattering generation the 
process would be repeated (with slightly different energies, angles etc.). Consequently, 
counter- streaming is maintained - at least for the duration of the prompt burst; the CBM 
slowly 'eats' momentum anisotropy from the GRB pulse. 

Whether the magnetic fields generated by the passing GRB photon burst are quasi- 
static to the extent that a trailing shocked ejecta could be influenced by the plasma and 
field structures produced, has yet to be determined. Field structures are sustained, we 
speculate, at least as long as significant photon momentum anisotropy is present. 



5. DISCUSSION AND CONCLUSIONS 



In conclusion, we have presented a new split dynamics framework - the Photon- 
PLASMA code - for modeling the plasma physics of hot tenuous plasmas. Large scale 
fields are modeled by means of a PIC code. Collisions and other sub-Debye processes, 
are modeled through a detailed Monte-Carlo approach. 

Our first wakefield simulation results using the code suggest that magnetic fields and 
plasma filaments are created in the wakefield of a prompt GRB. Particles are heated and 
accelerated and photons scatter into a distribution capable of photon-photon production. 

We identify two questions to be answered in the future, before significant progress 
can be made on GRB-CBM modeling and relativistic afterglow shock modeling: 

1. Can the plasma in the GRB wakefield undergo magnetic hysteresis - i.e. can it 
become magnetized, and stay magnetized, on timescales comparable to the time 
gap between the afterglow shock front and the GRB wakefield? 

2. Which binary particle processes, other than Compton scattering, are inclusions 
imperative to more further investigation of the effects found in the present work? 

The first of these questions is likely to require the extension of the modeling frame- 
work to include the transition from Vlasov-kinetic (PIC) codes to a coarser approxima- 
tion in the plasma physics formulation hierarchy (i.e. Vlasov-ffuids, Hall-MHD, . . . ). 

Regarding the second question; from figure SI right panel, we see that counter- 
streaming pairs of photons (above and below the green lines) are above pair production 
threshold. Such photon pairs are not observed in similar runs without the plasma effects 
included. We therefore argue that long-time development of the GRB-CBM system must 
include both plasma effects and photon-photon pair production. 

The energetic photons in the spikes of fig.|4](right panel) will, if they survive travers- 
ing the CBM, be observed as non-thermal emission. Although the color of these spikes 
indicate weights only of order 10~ 4 and less, it should be noted that these spiky struc- 
tures accumulate with time. They may eventually after passing a much longer distance 
lead to significant influence on the GRB spectrum away from its initial profile. 

Even if the interaction seemingly is locally unimportant for the burst photon front, 
such plasma and pair production effects could eventually affect - in an accumulative 
way - the final resulting observed GRB prompt burst. 
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